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Within the scope of an anisotropic Bianchi type-VI cosmological model we have studied 
the evolution of the universe filled with perfect fluid and dark energy. To get the determin- 
istic model of Universe, we assume that the shear scalar (a) in the model is proportional 
to expansion scalar (#). This assumption allows only isotropic distribution of fluid. Exact 
solution to the corresponding equations are obtained. The EoS parameter for dark energy as 
well as deceleration parameter is found to be the time varying functions. Using the observa- 
tional data qualitative picture of the evolution of the universe corresponding to different of 
its stages is given. The stability of the solutions obtained is also studied. 
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I. INTRODUCTION 

After the discovery of late time accelerating mode of expansion of the Universe a number of 
models are offered to explain this phenomenon. Most of the dark energy models such as cosmo- 
logical constant, quintessence, Chaplygin gas etc. are modeled with a constant EoS parameter. 
Recently in a number of papers different cosmological models with time dependent EoS param- 
eter was studied 01 |2l [TOJ [TT} [141 [20J. The aim of the current paper is to extend that study for a 
Bianchi type- VI cosmological model. 



H. BASIC EQUATIONS 

A Bianchi type-VI model describes an anisotropic but homogeneous Universe. This model was 
studied by several authors []3l[T2l[T3l[T5l[T6l[T9l, specially due to the existence of magnetic fields 
in galaxies which was proved by a number of astrophysical observations. 

Bianchi type-VI model given be given by fT2l [T3l 

d s 2 = dt 2 - a\e- 2mz dx 2 - a\e 2nz dy 2 - ajdz 2 , (2.1) 
with a\, a2, 6(3 being the functions of time only. H ere m, n are some arbitrary constants and the 



velocity of light is taken to be unity. The metric ( |2.1[ ) is known as Bianchi type-V I mo del. A 
suitable choice of m, n as well as the metric functions a\, ci2, a 3 in the BVI given by ( |2.1[ ) evokes 
Bianchi-type VIo, V, III, I and FRW universes. 

Here we consider the case when the energy momentum tensor has only non-trivial diagonal 
elements, i.e. 

rf=diag[r °,r 1 1 ,r 2 2 ,r 3 3 ] (2.2) 



Einstein field equations for the metric (2.1 ) on account of (2.2 ) have the form [[T2 
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We define the spatial volume of the model (2.1 ) as 



V = a\a2a-i, (2.4) 

and the average scale factor as 

a = V 1/3 = (a l a 2 a 3 y /3 . (2.5) 
Let us now find expansion and shear for BVI metric. The expansion is given by 

$ = u%=u»+T» a u a : (2.6) 
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and the shear is given by 



with 



where the projection vector P: 

P = P ■ PflV = §HV Uj^Uv-, Py = &V U^Uy. 

In comoving system we have = (1,0,0,0). In this case one finds 
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As one sees, neither the expansion nor the components of shear tensor depend on m or n., hence 
the Bianchi cosmological models of type VI, VIo, V, III and I has the same expansion and shear 
tensor. 

The Hubble constant of the model is defined by 



d 1 / d] d 2 d 3 ~ 
H=- = - (-+-+- 
a J \a\ a 2 a 3 . 
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The deceleration parameter q, and the average anisotropy parameter A m are defined by 
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where Hi are the directional Hubble constants 

d\ 
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(2.18) 



Note that, none of the above defined quantity depends on m or n, hence will be valid for not 
only BVI, but also for BVI , BV, Bill and BI. 
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III. SOLUTION TO THE FIELD EQUATIONS 

From ( |2.3e[ ) immediately follows 

fCl2\n , 

k\ = const. 



«3 «3 



(3.1) 



We also impose use the proportionality condition, widely used in literature. Demanding that 
the expansion is proportion to a component of the shear tensor, namely 



=jV 3 C7 3 3 . 



(3.2) 



The motivation behind assuming this condition is explained with reference to Thorne [18J, the 
observations of the velocity-red-shift relation for extragalactic sources suggest that Hubble expan- 
sion of the universe is isotropic today within w 30 per cent (6113. To put more precisely, red-shift 
studies place the limit 

^<0.3, (3.3) 

on the ratio of shear o to Hubble constant H in the neighborhood of our Galaxy today. Collins et 
al. (1980) have pointed out that for spatially homogeneous metric, the normal congruence to the 
homogeneous expansion satisfies that the condition § is constant. 
On account of (12. 101) and (12. 131) we find 



a 3 = N V 



A^o = const. 



(3.4) 



In view of (2.4) and (3.4) from (3.1 1 we find 
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Thus, we have derived metric function s in te rms of V. I n order to find the e quatio n for V we 
take t he following st eps. S ubtr actio ns of ( |2.3a[ ) from ( |2.3b[ ), ( |2.3c[ ) from ( |2.3c| ), and ( |2.3c[ ) from 
(|2.3a[) on account of (|3.5[), (|3.6|) and (|3.4[) give 
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where X12 = 3(m — n) /N^m + n), X23 = —3m/N^(m + n) and X31 = 3n/N?,{m + n). From ( |3.7[ ) 
immediately follows 
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After a little manipulation, it could be established that 
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Thus we conclude that under the proportionality condition, the energy-momentum distribution of 
the model should be strictly isotropic. 

Let us now gp back to the equation for V that now reads 



y-AoV^-6)/3* 3=0) Ao = ^(m + n 



which allows the solution in quadrature 

dV 
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Thus we have the solution to the corresponding equation in quadrature. 
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FIG. 1 . Evolution of the Universe given by a B VI cosmological model. 



Fig. |TJ shows the evolution of the Universe. As one sees, it is an expanding one. 



IV. PHYSICAL ASPECTS OF DARK ENERGY MODEL 

Let us no w find the e xpres sions for p hysical quantities. 

Inserting ( |3.11| ) into ( |2.15| ) and ( |2.16[ ) one finds the expression for expansion f}, Hubble param- 
eter H: 

(4.1) 
(4.2) 



t} = 3H= JAiV-^+V/^+Co/V 2 , 



and deceleration parameter 
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The anisotropy parameter A m has the expression 
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Am= N 2 (m + n) 2 ' 



(4.3) 



6 



Bijan Saha 



The directional Hubble parameters are 
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FIG. 2. Evolution of the Hubble parameter 
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FIG. 3. Evolution of the deceleration parameter. 



Figs. and [|3j show the behavior of the Hubble parameter and deceleration parameter, re- 
spectively. 



From ( |2.3d| ) we find the expression for energy density 

e = T ° = - [x l V- 2 ~X 2 V-^+ 6 ^] , 
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(4.5) 
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where 



Further we obtain 
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FIG. 4. Evolution of the energy density 
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FIG. 5. Evolution of the EoS parameter. 
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Figs. Q and |[5j show the behavior of the energy density and EoS parameter, respectively. As 
we see, energy density is a decreasing function of time, while the EoS parameter changes its sign. 

So, if the present work is compared with experimental re sults obtained in |51 [8l |9l [TTJ, then 
one can conclude that the limit of CO provided by equation ( |4.6[ ) may accommodated with the 
acceptable range of EoS parameter. Also it is observed that for V = V C , 0) vanishes, where V c is a 
critical Volume given by 

< x \ 3N 3 /{4N 3 -6) 



Vr 



X 4 J 



(4.7) 



Thus, for this particular volume, our model represents a dusty universe. We also note that the 
earlier real matter at V < V c , where CO > later on at V > V c , where CO < converted to the dark 
energy dominated phase of universe. 

For the value of CO to be in consistent with observation [8], we have the following general 
condition 

Vi < V < V 2 , 



where 
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For this constrain, we obtain —1.67 < co < —0.62, which is in good agreement with the limit 
obtained from observational results coming from SNe la data [[8]|. 

For the value of CO to be in consistent with observation ifTTl . we have the following general 
condition 

V 3 <V<V 4l (4.11) 

where 



and 
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(4.12) 



(4.13) 



For this constrain, we obtain — 1.33 < CO < —0.79, which is in good agreement with the limit 
obtained from observational results coming from SNe la data [fTTl . 

For the value of CO to be in consistent with observation 0H1, we have the following general 
condition 

v 5 < V < v 6 , 



where 



and 



V* 



v 6 



Xi + 1.44Xi \ 



3Ay(4JV 3 -6) 



x 4 + 1.44X2 

X 3 + 0.92Xi \ 
X 4 + 0.92X2 



3^3/(4Af 3 -6) 



(4.14) 



(4.15) 



(4.16) 



For this constrain, we obtain —1.44 < co < —0.92, which is in good agreement with the limit 
obtained from observational results coming from SNe la data j51|9l- 
We also observed that if 

2Xi \ 3^3/(4^3-6) 

Vq 



X 4 +X2 



(4.17) 
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then for V = Vq we have ft) = — 1, i.e., we have universe with cosmological constant. If V < Vq 
the we have CO > — 1 that corresponds to quintessence, while for V > Vo we have CO > — 1, i.e., 
Universe with phantom matter |@J. 

From ( |4.5[ ) we found that the energy density is a decreasing function of time and £ > when 

,<(§) . (4,8, 

In absence of any curvature, matter energy density Q. m and dark energy density Q.^ are related 
by the equation 

Q m + a A =-^ + — -^ = 1. (4.19) 



3H 2 3H 2 



Inserting ( |4.1[ ) and ( |4.5[ ) into ( |4.19[ ) we find the cosmological constant as 

A = [3Cl - (Xi/x^V- 2 + [3Ai -X 2 /k]V- 2{n ' +3 ^ 3n \ 



(4.20) 



As we see, the cosmological function is a decreasing function of time and it is always positive 
when 

'X l /K-3C \ 3N ^- 6) 



V > 



3Ai-X 2 /k) 



(4.21) 



Recent cosmological observations suggest the existence of a positive cosmological constant A 
with the magnitude A{Gfi/c i ) ~ 10~ 123 . These observations on magnitude and red-shift of type la 
supernova suggest that our universe may be an accelerating one with induced cosmological density 
through the cosmological A-term. Thus, the nature of A in our derived DE model is supported by 
recent observations. Fig. Q shows the evolution of the cosmological constant. As is seen, it is a 
decreasing function of time. 
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FIG. 6. Evolution of the cosmological constant 



For the stability of corresponding solutions, we should check that our models are physically 
acceptable. For this, the velocity of sound is less than that of light, i.e., 

0<v, = ^<1. (4.22) 
de 

In this case we find 
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dp _ X { - [(N 3 + 3)X 4 /3^V 3 ]y( 4jV 3-6)/3iV3 
Vs ~de~ Xi - [(N 3 + 3)X 2 /3N 3 ] y(4JV 3 -6)/3JV 3 
Fig. Ijvj shows the behavior of v s in time. 



(4.23) 
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FIG. 7. Speed of sound with respect to cosmic time 



As one sees, there are regions, where the solution is stable. Choosing the problem parameters, 
such as m, n,N 3 we can obtain the stable solutions. 



V. CONCLUSION 

In this report we have studied the evolution of the universe filled with dark energy within the 
scope of a Bianchi type- VI model. In case of a BVI model we found the exact solutions to the field 
equations in quadrature. It was found that if the proportionality condition is used, this together with 
the non-diagonal Einstein equation leads to the isotropic distribution of energy momentum tensor, 
i.e., Tj 1 = T 2 2 = T 3 3 . This fact allows one to solve the equation for volume scale V exactly. The 
behavior of EoS parameter CO is thoroughly studied. It is found that the solution becomes stable as 
the Universe expands. 
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